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Abstract 

In this letter we analyse an amazing variational approach to chemical reactions. Our 
results clearly show that the variational expressions are unsuitable for the analysis 
of empirical data obtained from chemical reactions. 

Some time ago Liu and He [1] proposed the application of a variational ap- 
proach to chemical reactions. Their method consisted of the transformation of 
the differential equation for the reaction rate into a kind of Newton equation 
that they easily converted into a variational problem for the mechanical en- 
ergy. By means of an appropriate ansatz they derived approximate expressions 
for the extent of reaction and half-time. 

The purpose of this letter is to investigate if those variational expressions may 
be of any utility for a chemist. 

Liu and He [1] considered the chemical reaction 
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If NA(t), Nsit), and Nc(t) are the number of molecules of the species A, B, 
and C, respectively, at time t then Liu and He [1] assumed that Na{0) = a, 
and Nb{0) = Nc{0) = 0. If we call x = Neit) = Nc{t), then we conclude that 
NA{t) = a — nx, where x is known as the extent of reaction [2]. The unique 
rate of reaction can be defined in terms of the extent of reaction as f = dx/dt. 
Liu and He [1] further assumed that the rate law is given by 



At this point we stress the fact that this expression is correct only if the chemi- 
cal reaction ([1]) is elementary, otherwise the rate law may be more complicated. 
Most chemical reactions are not elementary and therefore the reaction order 
and molecularity do not necessarily agree, as discussed in any book on physical 
chemistry [2]or chemical kinetics [3]. What is more, the order of reaction may 
not even be a positive integer [2,3]. This fact is known by any undergraduate 
student of chemical kinetics but Liu and He seem to be ignorant of it [1]. For 
concreteness here we assume that the rate law ([2]) is correct. 

Liu and He [1] tried and solved the differential equation ([2]) approximately by 
means of a variational approach based on the so-called semi-inverse method. 
It consisted in finding the minimum of the integral expression 



where 77 is a variational parameter. Notice that x^"*" (t) satisfies the boundary 
conditions at t = and t 00. They found that the optimal value of the 
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by means of the variational ansatz 
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effective first-order rate constant rj was given by 

V = (5) 



Liu and He [1] argued tliat cliemists and teclinologists always want to know tlie 
lialf-time ti/2 = t{x = a/2) (wliicli tliey decided to call halfway). According 
to equations (j4]) and ([5]) the half-time is given approximately by 



which they showed to be 98% accurate for n = 2. 

Any textbook on physical chemistry [2] or chemical kinetics [3] shows that the 
exact solution to equation ([2]) is 

1 [l + A;(n-l)a«-it]'/("-')J 
and that the exact half-time is given by 

j.exact fQ\ 

~k{n-l)a--^ 

Any undergraduate student of chemical kinetics is able to obtain the exact 
solution ([7j) to the differential equation ([2]) and, consequently, the exact half- 
time ([8]). Therefore, there is no need for approximate methods for this ex- 
tremely simple problem. However, we may assume that Liu and He [1] chose 
this exactly solvable problem as a benchmark for testing their powerful varia- 
tional approach. We will therefore analyse if their approximate result is worth 
the effort. 

It is common practice in chemistry to estimate the half-time from experi- 
mental data in order to determine the order of the reaction. Obviously, an 
inaccurate expression would lead to an inexact order of reaction. 



3 



The variational half-time ([6]) is reasonably accurate for n = 2 because it is 
exact for n = 1. The reason is that the variational ansatz (jlj) is the exact 
solution for a first-order reaction when i] = k. Notice that it leads to such a 

var /j-exa 
1/2/ ''1/2 

unity as n increases. 



result when n = 1. Table 1 shows that t^%/^T/T^ increasingly deviates from 



The half-time (or half-life) is a particular case of partial reaction times. We 
may, for example, calculate the time t = ti/4 that has to elapse for the number 
of A molecules to reduce to a/4 {x = 3a/4). It is not difficult to verify that 

^exact 

= 2"-^ + 1 (9) 

4-exact ^ ' 

From the experimental measure of ti/2 and ti/4 chemists are able to obtain 
the reaction order n. However, if they used the variational expression (jlj) they 
would obtain 

^ = 2 (10) 

4-var ^ ' 

that is useless for n ^ 1. According to what we have just said it is not surpris- 
ing that this ratio is exact for n = 1. We clearly appreciate that the variational 
result does not provide the information that chemists would like to have be- 
cause it only predicts first-order reactions. Most probably Liu and He [1] were 
ignorant of this fact when they proposed their remarkable variational method. 

From the discussion above we conclude that no chemist will resort to Liu and 
He's variational expressions in the study of chemical reactions. There is no 
reason whatsoever for the use of an unreliable approximate expression when 
one has simple exact analytical expressions at hand. Besides, we have clearly 
proved that the variational expressions are utterly misleading. 

In a review on asymptotic methods for strongly nonlinear problems He [4] 
came back to this problem and stated that the exact reaction extent for n = 2 
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x^^fin = 2) = a (l 
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This result is obviously wrong because it exhibits an unphysical pole at t = 
l/{ka). From this incorrect expression He [4] derived a meaningless negative 
half-time 



In order to obtain a reasonable agreement with the variational result ([6]) He [4] 
carried out the following wrong calculation 



In this way He [4] managed to obtain two unphysical negative half-times that 
agreed 98%. 

If you think that a paper with such a poor scientific quality is uncommon in 
the physics literature I will prove you wrong. There has recently been great 
interest in pseudo-scientific nonsensical results as we have discussed in several 
communications [5-13]. The journals that publish such extremely low quality 
papers do not accept comments or criticisms to them. I recommend the reader 
to have a look at an approach to population dynamics that predicts a negative 
number of rabbits [9]. There must be considerable gains in publishing such 
papers that compensate the obvious loss of credibility. 
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Table 1 

Ratio ^1/2/^1/2'^* in terms of the reaction order n 



n 


fvar uexact 
''l/2/''l/2 


2 


0.98 


3 


0.80 


4 


0.59 


5 


0.41 


6 


0.27 


7 


0.17 


8 


0.11 


9 


0.07 


10 


0.04 
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